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> 1 INTRODUCTION 



ABSTRACT 

It is expected that an average protostar will undergo at least one impulsive in- 
teraction with a neighbouring protostar whilst a large fraction of its mass is still in a 
massive, extended disc. If protostars are formed individually within a cluster before 
falling together and interacting, there should be no preferred orientation for such in- 
teractions. As star formation within clusters is believed to be coeval, it is probable 
that during interactions, both protostars possess massive, extended discs. 

We have used an SPH code to carry out a series of simulations of non-colpanar 
disc-disc interactions. We find that non-coplanar interactions trigger gravitational in- 
stabilities in the discs, which may then fragment to form new companions to the 
existing stars. (This is different from coplanar interactions, in which most of the new 
companion stars form after material in the discs has been swept up into a shock layer, 
and this then fragments.) The original stars may also capture each other, leading to 
the formation of a small- A/" cluster. If every star undergoes a randomly oriented disc- 
disc interaction, then the outcome will be the birth of many new stars. Approximately 
two-thirds of the stars will end up in multiple systems. 
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Interactions between young protostars are believed to play 
an important role in the star formation process, with partic- 
ular implications for the forma tion o f binary systems, and for 
the lifetimes of discs. Larson (199C) suggested that capture 



during interactions might be an important binary formation 
mechanism, but further investigations (Clarke & Pringle 



1991a, jJlarke & Pringle 1991b, Hall, Clarke & Pringle 1995, 
Heller iyyo) indicated that capture was unlikely to be a 



rounded by massive, extended discs, and that most of these 
interactions may be prograde and coplanar. If, however, star 
formation is more distributed, and interactions occur when 
two protostars formed separately within a cluster fall to- 
gether, then there should be no preferred orientation for 
interactions. 

This is the last of three papers describing th e results 
of simulation s of protostellar interactions. Paper I (Walking 



dominaint meclianism. 



et al. 1997c) described interactions in whi ch only one of 
the st ars possessed a disc, and Paper II (Watkins et al 



In all of the above cases, only interactions in which just 
one of the stars possessed a disc were investigated. How- 
ever, observations indicate that stars within clusters tend to 
be formed coevally (e.g. Gauvin & Strom 1992), and so it 
might be expected that both interacting stars would possess 
extended discs. In particular, if star formation is triggered 
dynamically, then the stars formed are likely to have a small 
spread of ages, and so again disc-disc interactions should 
tend to predominate over star-disc interactions. 

Numerical simulations of star formation triggered by 
the collision of two clumps within a molecular cloud (Turner 
et al. 1995, Whitworth et al. 1995) suggest that the major- 
ity of protostars will undergo interactions whilst still sur- 



1997b) presented the results of simulations of coplanar in- 



teractions in which both protostars had massive, extended 
discs. This paper investigates non-coplanar interactions be- 
tween two stars each possessing a massive, extended disc, 
analyses the results presented in both this paper and Paper 
II, and compares them with the results for star-disc interac- 
tions presented in Paper I. 



2 PHYSICAL AND COMPUTATIONAL 
MODEL 

The numerical method used is described in Paper I. The 
discs modelled in the simulations described in this paper 
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property 



value 



disc radius 

disc mass 

star mass 

disc density profile 

eccentricity of orbit 

initial separation 

shear viscosity 



lOOOAU 

O.5M0 

O.SMq 

,.-3/2 

1.0 

5000AU 

Vs = Csh/800 



Table 1. Physical parameters used in simulations 



Table 2. 


List of simulations 




Run 


^'peri 


4> 


dd05-08 


500, 1000, 1500, 2000 


(f) = 7r/4 


dd09-12 


500, 1000, 1500, 2000 


<j> = tt/2 


ddl3-16 


500, 1000, 1500, 2000 


4> = 3tt/4 



have the same physical properties as the discs used in Paper 
I. A brief summary of these properties is given in Table 
^ Each disc is modelled with 2000 SPH particles, so it is 
not resolved vertically. One particle represents a mass of 
2.5 X 10~*Mq. 

Non-coplanar disc-disc interactions are simulated at the 
same periastra and orbit inclinations as the star-disc encoun- 
ters, in order to limit the parameter space to be searched. 
In all cases, the spin of the secondary disc is parallel to the 
orbit. A total of 12 simulations of non-coplanar disc-disc 
interactions has been carried out. 

The simulations are labelled so that the number of the 
disc-disc simulation is the same as that of the corresponding 
disc-star simulation from Paper I, so that for instance the 
disc-disc simulation dd05 has the same orbital inclination 
and periastron as the disc-star simulation ds05. 



3 SIMULATIONS 

The figures presented in this section show the results of three 
of the simulations conducted. The figures are grey-scale plots 
of surface density, in which the shading is logarithmically- 
scaled. These plots are overlaid with contours of constant 
surface density, that are equally separated in log-space. Oc- 
casionally an extra contour has been added in order to high- 
light an important feature of the figure, such as a shock or 
spiral arm. The contour levels used are given in the caption 
to the figure. The time at each figure, in units of lO** years, 
is also given, and is calculated from the start of the simula- 
tion, when the two stars are at a separation of bvdisc- The 
individual frames within a figure, when referred to within 
the text or a caption, are labelled from top left to bottom 
right in rows, so that frames 1 and 2 are on the top row, 
frames 3 and 4 on the next row and so on. 



of the primary, and passes up through it. There is no sig- 
nificant shock as the secondary disc penetrates through the 
primary. Instead the evolution of the system is dominated 
by a tidal-arm instability. A spiral arm forms joining the 
secondary disc to the primary star, and likewise joining the 
primary disc to the secondary star. Corresponding opposing 
tails form in the discs. The arms joining the stars are quickly 
dissipated and merge with the central discs surrounding each 
star, while the tidal tails remain in a circumbinary disc-like 
distribution. The inner regions of the discs undergo strong 
spiral-arm instabilities that lead to their fragmentation. The 
disc surrounding the primary produces one new condensa- 
tion, whilst that around the secondary produces three. Two 
of these fragments merge, so that at the end of the simula- 
tion the primary is in a binary and the secondary is in an 
unstable triple system. The primary has a circumstellar disc 
of mass O.18M0 and radius 160AU, whilst its companion 
is of mass 0.06Mq. The binary orbit has eccentricity 0.14 
and periastron 185AU. The secondary has a circumstellar 
disc of mass 0.24Mq and radius llOAU. The companions 
are accreting from this disc and truncating it. The compan- 
ions have masses 0.04Mq and O.OOMq. The binary and the 
triple system are bound to each other, with an eccentricity 
of 0.7 and a periastron of 860AU. The primary disc has been 
twisted through an angle of 9° , and its binary companion is 
orbiting in the plane of the disc. The secondary disc has 
been twisted though 5°, and likewise its companions lie in 
the plane of the disc. 



3.2 run ddlO : </> = 7r/2 

Figure ^ shows the results for run ddlO, in which the perias- 
tron is lOOOAU. The secondary is initially below the primary, 
and orthogonal to it. As the secondary slices up through the 
primary, the leading edge of the secondary becomes decel- 
erated. This triggers fragmentation of the secondary disc to 
form two new condensations. These merge to form a binary 
companion to the secondary star. As this binary orbits, it 
sweeps out a gap in the surrounding material. In the mean- 
time, the interaction has caused the primary disc to frag- 
ment to produce four new condensations. Figure ^ shows 
the fragmentation of the primary disc viewed from above. 
The two frames correspond to the last two frames of Fig. ^ 
The secondary has a circumstellar disc of radius 30AU 
and mass 0.16Mq. Its binary companion is of mass 0.14Mq, 
and their orbit has eccentricity 0.16 and periastron 90AU. 
The primary has a tidally-truncated circumstellar disc of 
mass O.ISMq and radius 60AU, and is in a small- A/" cluster 
with companions of mass 0.02, 0.04, 0.06 and 0.07Mq. The 
cluster and the binary are bound to each other, with eccen- 
tricity 0.19 and periastron 930AU. The circumstellar disc 
around the primary is untilted by the encounter, but the 
outer disc and the orbits of the newly-formed companions 
have been twisted by approximately 12°. 



3.1 run dd06 : (j) = 7r/4 

In this simulation, the orbit and the spin of the secondary 
are inclined at 45 degrees to the spin of the primary, and the 
periastron is lOOOAU. The secondary starts below the plane 



3.3 run ddl3 : (p = 37r/4 

In runs ddl3-16, the orbit and the secondary disc are in- 
clined at 135 degrees to the primary disc. The results of run 
ddl3 are shown in Fig. The periastron is 500AU. The 
secondary is initially below the primary. As it penetrates 
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Figure 1. Simulation dd06 : 4> = 'r/^, fperi = lOOOAU. 2000 x 1500 AU region, x — y plane. Contour levels at [0.3, 3, 9 (frame 3, 4), 
30, 300] g cm~-^. The primary is initially on the left of the figure, rotating clockwise and in the plane z = 0. The secondary is initially 
below and to the right of the primary. Both the orbit and the secondary are inclined at 7r/4 to the primary. Viewed from this angle, the 
orbit is in a clockwise sense. As they interact, spiral arm instabilities cause the fragmentation of the discs, resulting in the formation of 
a binary and a triple system. 
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Figure 2. Simulation ddlO : cf) = 7r/2, Vperi = lOOOAU. 2000 x 1500 AU region, y - z plane. Contour levels at [0.3, 2 (frame 6), 3, 
6 (frame 2), 7 (frame 3, 4), 30, 300] g cm~^. The primary is initially at the top of the figure and is viewed edge-on. The secondary is 
initially below the primary. Viewed from this angle, the secondary is rotating clockwise, and the orbit of the encounter is also clockwise. 
The secondary disc fragments to form a binary, whilst the primary disc fragments to produce four new low-mass protostars. 
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Figure 3. Simulation ddlO : 4> = 7r/2, Vp^ri = lOOOAU. 1000 X 750 AU region, x - y plane. Contour levels at [0.3, 3, 30, 300] g cm"^. 
View from above showing fragmentation of the primary disc to produce four new protostars. 



Figure 4. Simulation ddlO : (j> = n/2, Vp 



lOOOAU. Positions of SPH particles at end of simulation. The left-hand figure shows the 



primary star and its four companions. The right-hand figure shows the secondary star, its binary companion and the circumbinary disc 
in which they have cleared a gap. 



through the primary disc, spiral arms form which join the 
secondary disc and the primary star, and likewise the pri- 
mary disc and secondary star (frame 2). Due to the relative 
inclinations of the discs, these arms both lie on the same 
side. The outer regions of the discs then form into trailing 
arms, while the original arms dissipate as the stars move 
apart (frame 3). The central discs around each star undergo 
strong spiral-arm instabilities without fragmenting (frames 
4 & 5), before settling down to a steady state (frame 6). At 
the end of the encounter, the primary has a circumstellar 
disc of mass 0.44Mq and radius 60AU, and the secondary 
has a disc of mass O.STMq and radius 125AU. The two stars 
are in a binary orbit with eccentricity 0.32 and periastron 
400 AU. 



4 RESULTS 

In coplanar disc-disc encounters, disc material is swept up 
into a shock layer that then fragments to produce new pro- 
tostars. The remaining discs may later undergo more frag- 
mentation via spiral arm instabilities triggered by the en- 
counter. In the case of non-coplanar encounters, shocks do 
not play a major part in the evolution of the system. Instead, 
the spiral-arm instabilities caused by the interaction are the 
sole source of fragmentation and energy transfer from the 
orbit of the discs. 

Frames 1-3 of Fig. ^ show the mass that ends up bound 
to each of the original stars, and the mass that becomes un- 
bound from the system, for each of the non-coplanar encoun- 
ters. The mass that is bound to a star may be either in the 
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Figure 5. Simulation ddl3 : </> = 37r/4, rperi = 500AU. 2000 x 1500 AU region. Contour levels at [0.3, 3, 9 (frame 2, 4), 30, 300] g cm-^. 
The primary is initially at the right-hand side of the figure, rotating clockwise and in the plane z = 0. The secondary is initially below 
and to the left of the primary. Both the orbit and the secondary are inclined at 37r/4 to the primary. Viewed from this angle, the orbit is 
in an anti-clockwise sense. The encounter does not lead to fragmentation of the discs, but the primary and secondary are captured into 
a binary system. 
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4> = -n/4: 



bound to primary 
boLiiici Lo secondary 
unbound 



1.5 



0.5 



= ^/2 



bound to primary 
bound to secondary 
unbound 



0.5 



1.5 



0.5 



1 1.5 

rp.ri / I'd,,, 



0=3^/4 
' I ' ' ' ' I 
bound to primary 
bound to secondary 
unbound 




-10 



0.5 



-15 



0.5 



Figure 6. Frames 1-3 show the fraction of the disc mass that ends up bound to the primary star, bound to the secondary star and 
unbound, for (j> = 7r/4, it/2 and 37r/4 encounters respectively. Frame four shows the change in energy of the original orbit for all 
non-coplanar encounters. 



form of a circumstellar disc, or bound up in newly-formed 
companions. The results for the corresponding disc-star en- 
counters are given in Paper I. A similar amount of mass is 
unbound from the system for the star-disc and disc-disc en- 
counters, but in the disc-disc encounters the secondary does 
not capture any material from the primary. 



4.1 Energy exchange 

Frame four of Fig. ^ shows the change in the energy of the 
orbit of the discs due to the encounter. The energy change is 
given in code units, in which the total binding energy of the 
disc is 6.6 and the binding energy outside 0.5rdisc is 0.82. 
As expected, the encounters with </> = 37r/4, which have a 



retrograde element to the orbit, remove more energy from 
the orbit than those with — 7r/4. However, the most dis- 
sipative non-coplanar encounters, as well as being the ones 
that cause the most fragmentation, are the orthogonal en- 
counters. This is in contrast to the corresponding disc-star 
encounters, in which the orthogonal encounters remove the 
least energy from the orbit. 

The energy removed from the orbit is much greater for 
disc-disc encounters than for the corresponding star-disc en- 
counters, the results for which are shown in Fig. 6 of Paper 
I. For encounters with </!> = 7r/4, the disc-disc encounters 
remove 2-3 times as much energy from the orbit as do the 
corresponding disc-star encounters; for those with = 7r/2 
and (j> = 37r/4, the ratio reaches as high as 20. Due to the 
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- 
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Figure 7. Distribution of masses of newly-formed stars, and orbital eccentricities as a function of periastron after disc-disc encounters 



alignment of the spin of the secondary disc with that of 
the orbit in all of the encounters presented here, energy will 
be transferred to the orbit via a corotation resonance. It is 
therefore likely that the results here provide a lower limit to 
the energy dissipated from the orbit during such encounters. 
Encounters with random inclinations, in which there is no 
resonance, will remove even more energy from the orbit of 
the stars. As a result of this, capture is likely to be a more ef- 
ficient binary formation mechanism for disc-disc encounters 
than for star-disc encounters. 

4.2 Properties of new stars 

As well as removing energy from the orbit of the stars, the 
gravitational instabilities in the disc triggered by the en- 
counter can lead to fragmentation of the disc and the for- 
mation of new companions to the original stars. For the 
non-coplanar star-disc encounters, fragmentation occurred 
in just one out of twelve simulations, leading to the forma- 
tion of two companions to the primary star. For the corre- 
sponding disc-disc encounters, fragmentation occurs in eight 
of the 12 cases, producing a total of 16 new stars. Of those 
16 stars, seven are produced by fragmentation of the pri- 
mary disc, and nine by fragmentation of the secondary, so 
that it seems there is no particular preference for which disc 
fragments. 

The left-hand frame of Fig. ^ shows the final distribu- 
tion of the masses of the protostars that are formed during 
the disc-disc encounters. Most of the fragments form with 
masses of around 0.04Mq. Condensations that form by frag- 
mentation of the shock layer in coplanar encounters accrete 
from the discs and end up with masses of approximately 
0.12Mq, whilst the condensations formed in the SOAQ en- 

* SOA stands for spin orbit antiparallel, i.e. the two spins are 
parallel to one another and antiparallel to the orbit (see Paper 
II). 





SOP 


^=1 




<^= ^ 


SOM 


SOA 


500 


2 





1 





3 





1000 


3 


3 


5 


1 


7 


2 


1500 


3 





2 


1 


3 


2 


2000 


2 


1 


2 





1 


2 



Table 4. Number of companions formed by disc fragmentation 

counters accrete the most material and end up with masses 
typically of 0.2Mq. In most cases there is still disc material 
that will be accreted, and in a physical situation there will 
also be a remnant infalling envelope, so that this distribu- 
tion of masses may evolve greatly, and therefore may not 
be particularly close to the final outcome. The right-hand 
frame of Fig. M shows the distribution of eccentricities as a 
function of periastron for all of the systems that are formed 
during disc-disc encounters. It can be seen that the encoun- 
ters tend to change the eccentricities but not the periastra 
of the original orbits, leading to a clustering of points at 
rperi/rdisc=0-5, 1.0, 1.5, 2.0. The orbits of the stars that 
are formed during the encounter have smaller eccentricities 
and periastra. Although there is a lot of scatter in the data, 
a general trend can be seen from having low eccentricities 
at low periastra to high eccentricities at larger periastra, 
as is observed amongst both pre-main-sequence and main 
sequence stars (Mathieu 1994, Duquennoy & Mayor 1991). 

4.3 Capture rates 

In order to attempt to ascertain how important capture and 
fragmentation are for the formation of multiple systems, the 
calculations of Paper I are repeated for the disc-disc en- 
counters. We calculate both the capture rate and the rate 
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Region 


n/pc ^ 


y/km 


Tc/Myr-i 


rl 


/Myr-l 


Trapezium (centre) 


10^ 


1.5 


0.6 




0.03 


Trapezium 


2 X 10^ 


1.5 


0.12 


6 


X 10-3 


Open cluster 


10^ 


1 


0.013 


2 


X 10-^ 


GMC 


12 


2 


3 X 10-4 


5 


X 10-6 



Table 3. Capture rates for disc-disc encounters in different star-forming environments 




V / km s-i 

Figure 8. Capture rate per unit number density for disc-disc 
encounters, as a function of velocity dispersion. The solid line 
shows the rate for randomly distributed inclinations, and the dots 
show the rate corrected for previous disruption of discs. If all 
encounters are SOP then the capture rate is negligible. 




0.5 1 15 a 2.5 

V / km s-i 



Figure 9. Mean probability of capture during a disc-disc en- 
counter for a randomly-inclined encounter, as a function of veloc- 
ity dispersion. 



of formation of companions by fragmentation. As fragmen- 
tation occurs for almost all of the disc-disc encounters, the 
rate of capture events will not in general be the rate of for- 
mation of binaries by this mechanism, but rather the rate 
of capture of single stars and multiple systems into hierar- 
chical configurations. We consider only encounters in which 
the spin of the secondary disc is parallel to that of the orbit, 
so that for a coplanar encounter with inclination tf) — n, we 
use the SOmQ rather than the SOA results. 

Figure ^ shows the capture rate per unit number density 
for disc-disc interactions, Fc, as a function of velocity disper- 
sion. It filao ohowo the rate corrected for proviouo oncountoro. 



Fc whidh assumes that an encounter with r^^, j < 2raij^ com- 



pletely destroys both discs. This information is then used 
to calculate capture rates in the star-f orming environments 
considered in Clarke & Pringle (1991a) and in Paper L The 
rates are given in Table ^ It can be seen that in typical star- 



t SOM stands for spin orbit mixed, i.e. one disc spins parallel to 
the orbit, and the other spins antiparallel to the orbit (see Paper 
II). 



forming environments such as an open cluster or a GMC, 
capture does not happen at a significant rate, mainly due 
to the low number of encounters that occur. The rates cor- 
rected for previous encounters are also much lower than the 
uncorrected rates, due to the fact that disc-disc encounters 
destroy the discs out to a greater periastron than disc-star 
encounters. 

If star formation is dynamically-triggered, however, for 
instance by the collision between two clumps within a molec- 
ular cloud, then it is expected that almost every protostar 



form ed will undergo at least one interaction (Turner et al 



1995). The probability of capture during an encounter is 
then pc — Fc/Fh, where F^ is the rate of encounters with 
rperi < 2rdisc- Figure ^ shows Pc as a function of the velocity 
dispersion. In approximately 10% of cases, the stars will be 
bound to each other after the encounter. Given that frag- 
mentation occurs in most of the disc-disc encounters, it will 
usually not be two individual stars that become bound to 
each other in a binary following capture. Rather, at least one 
of the two stars will already have one or more companions, 
so that the outcome will be a hierarchical system. 
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Region 


n/pc ^ 




r//Myr-l 


r^/Myr-l 


Tp/Myr-i 


Trapezium (centre) 


10* 


1.5 


3.0 


2.0 


4.8 


Trapezium 


2 X 10^ 


1.5 


0.6 


0.4 


1.0 


Open cluster 


10^ 


1 


0.05 


0.03 


0.08 


CMC 


12 


2 


2 X 10"^ 


1 X 10"^ 


3 X 10^3 



Table 5. Rate of formation of companions for disc-disc encounters in different star-forming environments. is the uncorrected rate 
assuming random inclinations for encounters, V ^ is the corresponding corrected rate, and Fp is the uncorrected rate if all encounters are 
coplanar prograde. 




1.5 



0,5 - 



0.5 



1.5 

km a" 



2.5 



Figure 10. Rate of formation of companion stars, per unit num- 
ber density for disc-disc encounters, as a function of velocity dis- 
persion. The solid line shows the rate for randomly distributed 
inclinations, the long dashes show the rate assuming all encoun- 
ters are prograde coplanar, and the short dashes show the rate 
corrected for previous disruption of discs. 



If all of the interactions are sop| encounters, then the 
capture rate is negligible, due to the resonances that tend 
to unbind the stars from each other. 



4.4 Fragmentation rates 

The rate of formation of companions has also been calcu- 
lated using the same method as in Paper I. Table ^ lists 
the number of companions formed during each of the disc- 
disc encounters simulated. When calculating the rate of for- 
mation of companions, we divide these numbers by two to 
obtain the number of companions formed per disc in the 
simulations. Figure ^ shows the rate of formation of com- 
panions per unit number density. This information is then 



■t SOP stands for spin orbit parallel, 
to the orbit (see Paper II). 



i.e. both discs spin parallel 



Figure 11. Mean number of companions formed in a disc-disc en- 
counter, as a function of velocity dispersion. The solid line shows 
the number for randomly distributed inclinations, the dots show 
the number assuming all encounters are prograde and coplanar. 



used to calculate the companion formation rates in typi- 
cal star- forming regions, given in Table |^. In high-density 
regions like the Trapezium, fragmentation during disc-disc 
collisions can produce a high number of companions, so that 
for instance in the very core of the Trapezium, a star will 
have on average 2 to 5 companions after IMyr. In lower den- 
sity regions, the rates are much lower due to the low number 
of interactions. 



5 DISCUSSION 

5.1 Disc-star interactions 

Disc-star interactions provide a mechanism for removing en- 
ergy from the orbit of the two interacting stars, for truncat- 
ing the disc, and for triggering the fragmentation of the disc 
to produce new stars. The results presented in Paper I are 
in general agreement with those of previous workers, with 
prograde coplanar interactions being the most destructive 
to the disc, but retrograde coplanar interactions removing 



© 0000 RAS, MNRAS 000, 000-000 



Non-coplanar disc-disc encounters 11 



the most energy from the orbit. The evolution of the pen- 
etrating, retrograde, coplanar interactions is dominated by 
the formation of a trailing shock behind the perturbing star 
that dissipates large amounts of energy and leads to the frag- 
mentation of the disc. For other penetrating interactions, 
the energy removed from the orbit is approximately twice 
the binding energy of the disc outside of periastron. 

If star formation within typical star-forming regions is 
distributed, then the capture of the two stars into a binary 
system due to the energy dissipated in an interaction hap- 
pens at a negligible rate, primarily due to the low number of 
discs that experience interactions. Even if every disc under- 
goes an interaction, as might happen during dynamically- 
triggered star formation, then the capture rate is extremely 
low unless the stars initially have a very low velocity disper- 
sion. If the majority of interactions are prograde, coplanar 
encounters, as predicted for stars that are formed from the 
same core or disc, or that are formed during the collision of 
two clumps within a molecular cloud, then the interactions 
will tend to unbind the stars, and no binaries will form by 
capture. 

Disc-star interactions, as well as removing energy from 
the orbit of the interacting stars, can also trigger the frag- 
mentation of the disc via gravitational spiral-arm instabili- 
ties, leading to the formation of one or two new stars that 
are bound to the original star. If every disc undergoes one 
disc-star interaction, then approximately 30% of stars will 
end up in multiple systems, rising to about 60% if all of the 
interactions are prograde and coplanar. Thus disc-star inter- 
actions can give rise to a large number of multiple systems, 
formed primarily by fragmentation of the disc triggered by 
the interaction, rather than by capture. 

5.2 Disc-disc interactions 

If the star formation within a region is coeval, as will be 
the case if it is dynamically triggered, then most encounters 
are likely to be disc-disc interactions, rather than star-disc 
interactions. For disc-disc interactions, the evolution of the 
system is almost always dominated by the fragmentation of 
the disc to produce multiple new companions to the origi- 
nal stars. This fragmentation occurs in two forms. The first, 
which happens for non-coplanar encounters, is via the same 
gravitational spiral-arm instabilities as for the star-disc en- 
counters. In the second mechanism, which dominates for the 
coplanar disc-disc encounters, the disc material between the 
two interacting stars is swept up into a shock layer that frag- 
ments to produce multiple new stars. Most often these form 
binary or triple systems with one of the original stars, but 
occasionally larger clusters can form. 

Disc-disc encounters form 2-3 times as many compan- 
ions as star-disc encounters. Capture will occur in approx- 
imately 15% of disc-disc encounters, a significantly higher 
rate than for star-disc encounters. Given the high frequency 
of fragmentation, in most cases capture will produce a small- 
N cluster of several stars, and so the process of disc-disc 
interaction can repeat itself hierarchically. 

5.3 Disc masses 

The calculations presented in this paper and its two compan- 
ions have demonstrated that under suitable circumstances 



(i) interactions between massive protostellar discs, and/or 
between a massive protostellar disc and a naked star, can 
spawn additional protostars from the disc material; and (ii) 
most of the resulting protostars are in binary or higher mul- 
tiple systems, with orbits smaller than the periastron of the 
original interacting pair. Therefore, if these processes can, 
starting with a relatively wide interacting pair, repeat them- 
selves in an hierarchical cascade, it may be possible to create 
ever closer systems, and hence to populate the entire ob- 
served range of binary separations. However, this can only 
happen if the protostellar discs involved remain massive, so 
that there is a sufficient reservoir of material from which to 
form new protostars. 

This requirement places at least two conditions on the 
evolution of protostellar discs. The first is that in an isolated 
protostellar disc the rate at which matter looses angular 
momentum and spirals onto the central star must be slow. 
Specifically, the time-scale tviac on which an isolated proto- 
stellar disc accretes onto its central star should be longer 
than the mean time tint between interactions. There seems 
no reason to suppose that this condition is not met routinely 
in star-formation regions, tvisc is essentially the time-scale on 
which angular momentum is redistributed in the disc, due 
to turbulent viscosity, or magnetic stresses, or global gravi- 
tational instabilities; tvisc is likely to be significantly longer 
than the spin period of the disc tapi„. By contrast, when clus- 
ters of protostellar discs form (by fragmentation) from the 
same collapsing core, ti„t should be comparable with igpin. 
Indeed, it is possible that most accretion onto the central 
star occurs as a result of non-linear instabilities triggered by 
interaction with another protostar; and that between inter- 
actions an isolated protostellar disc evolves secularly, with 
matter accreting only very slowly onto the central star, and 
the remaining disc becoming increasingly extended - and 
therefore increasingly prone to interaction - until it is in- 
volved in (another) impulsive interaction. 

The second condition is that protostellar discs must be 
replenished by infall from the surrounding gas. Otherwise 
the hierarchical cascade will tend to deliver protostars with 
ever-decreasing, and eventually sub-stellar, masses. Again 
this condition is likely to be fulfilled. Both theory and ob- 
servation seem to be agreed that the gravitational instability 
which creates the original protostellar discs is likely to be fol- 
lowed by continuing infall. As long as most of the infalling 
material has sufficient angular momentum relative to the in- 
dividual central stars, it will accrete onto their discs, thereby 
replenishing them in preparation for the next interaction. 

Thus, as long as (i) the evolution of an isolated proto- 
stellar disc is sufficiently slow, and (ii) protostellar discs are 
replenished by infall, the processes which we have modelled 
here could operate frequently, and in an hierarchical cascade 
to generate a wide range of binary (and higher multiple) sys- 
tems. 



6 CONCLUSIONS 

If star formation is not quiescent, but is dynamically trig- 
gered, then it is likely that almost every protostar undergoes 
at least one interaction whilst a large fraction of its mass is 
still in an extended disc. Figure hll shows the mean num- 
ber of companions formed in a disc-disc interaction. If every 
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protostellar disc undergoes an interaction which has a ran- 
dom inclination, then after the encounter the star wiU on 
average have 0.6 companions. If all of the interactions are 
coplanar SOP, then the mean number of companions is ap- 
proximately double that for randomly-inclined interactions. 
Some of the stars will be in binary systems, some in higher 
multiples, and some will stay as single stars. Many of the 
higher multiple systems will be in unstable configurations, 
and some stars will undergo slingshot events and be ejected 
from the birth site. 

For non-coplanar disc-disc interactions, the evolution 
of the system is dominated by spiral-arm gravitational in- 
stabilities. These instabilities can lead to the fragmentation 
of the discs, and the formation of new companions to the 
original stars. This occurs in two-thirds of the non-coplanar 
interactions modelled here. 

We have calculated the rates of capture due to disc- 
disc interactions. Given the large fraction of cases in which 
fragmentation occurs, any capture is likely to be of a mul- 
tiple system into a hierarchical grouping. We find that for 
randomly-oriented interactions, capture will occur in about 
20% of cases; if, however, all of the interactions are coplanar 
and spin-orbit parallel, then the rate of capture is negligible. 

We believe that the majority of protostellar discs will 
undergo at least one disc-disc interaction during their life- 
time. The results presented here indicate that such interac- 
tions lead to the formation of new protostellar discs in closer 
multiple systems. Therefore the compound effect could be 
to propagate the star-formation process from a few initial 
massive protostellar discs in wide binary systems, to lower 
masses and smaller separations, thereby populating most of 
the observed range of stellar masses and binary separations. 
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